Competing orders in the generalized Hund chain model at half-filling 
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By using a combination of several non-perturbative techniques - a one-dimensional field theoretical 
approach together with numerical simulations using density matrix renormalization group - we 
present an extensive study of the phase diagram of the generalized Hund model at half-filling. This 
model encloses the physics of various strongly correlated one-dimensional systems, such as two-leg 
electronic ladders, ultracold degenerate fermionic gases carrying a large hyperfine spin |, other 
cold gases like Ytterbium 171 or alkaline-earth condensates. A particular emphasis is laid on the 
possibility to enumerate and exhaust the eight possible Mott insulating phases by means of a duality 
approach. We exhibit a one-to-one correspondence between these phases and those of the two-leg 
electronic ladders with interchain hopping. Our results obtained from a weak coupling analysis are 
in remarkable quantitative agreement with our numerical results carried out at moderate coupling. 

PACS numbers: 71.10.Pm ; 71.10.Fd 



I. INTRODUCTION 

A major focus of the study of strongly correlated elec- 
tronic systems is the analysis of the competition between 
qualitatively distinct ground states and the associated 
quantum phase transitions (QPT)i in low dimensions. A 
reason to concentrate on these matters stems from the 
hope that the criticality of the system at such QPT's 
possibly results in a universal description of their vicin- 
ity. In 2 + 1 space-time dimensions, the resulting rela- 
tivistic quantum field theories which describe the zero- 
temperature transition between these quantum phases 
are in general strongly coupled and can be highly non- 
trivial^. 

In one dimension, the situation is much simpler since 
the quantum critical points in standard condensed mat- 
ter systems arc characterized by conformal field theories 
(CFT), which often admit a simple frcc-ficld represen- 
tation in terms of free bosons or fcrmions. In this re- 
spect, the bosonization approach has been very successful 
to investigate the physical properties of one-dimensional 
quantum phases^. Within this approach, several con- 
ventional and exotic long-range ordered phases have been 
revealed over the years in two-leg ladder models^"— and 
carbon nanotubc systemsi^Ti^^ at half-filling. Two dif- 
ferent classes of Mott-insulating phases have been found 
at half-filling in these systems. A first class is two-fold 
degenerate corresponding to the spontaneous breaking of 
discrete symmetries such as the translation symmetry in 
charge density wave (CDW) and bond-ordering phases, 
or the time-reversal symmetry in d-density wave (DDW) 
phase. In contrast, the second class of Mott-insulating 
phases is non-degenerate. A paradigmatic example of 
this class is the Haldane phase of the spin-1 Heisenberg 



chain2^ and of the two-leg spin ladder which breaks spon- 
taneously a non-local Z2 x Z2 discrete symmetry^ir— . 

Another striking particularity of many one- 
dimensional electronic systems is the existence of 
hidden duality symmetries, within the low-energy 
approach, which relate many of the competing orders to 
a conventional one like the CDW— li^i^"— . Actually, a 
general duality approach has been introduced recently 
to describe the zero-temperature spin-gapped phases 
of one-dimensional (ID) electronic systems away from 
half-filling2i. In this paper, we apply this approach 
to half-filled systems of four-component fermions and 
revisit the problem of competing orders in half-filled 
two- leg electronic ladders. In this particular case, some 
of the duality symmetries already exist at the level of 
the lattice model, and have been first revealed in Ref. [H. 
In addition to these, as it will be seen, there are also 
interesting emergent duality symmetries which relate 
non-degenerate Mott insulating phases to conventional 
order such as CDW. Unlike the former ones, those 
dualities do not bear a local representation on the 
lattice. 

The starting point of the duality approach to com- 
peting orders is to identify the internal global symmetry 
group H of the lattice model. For two-leg electronic lad- 
ders, or more generally two-band models, the building 
blocks of the model are four-component fermionic cre- 



ation operators on each site i: 



where / = 1 , 2 is the 



leg or orbital index and a =^, ^ denotes the spin-i index. 
Three basic global continuous symmetries are retained: 
a U(l) charge symmetry {cia.i e'''^cia,i), a SU(2) spin- 
rotational invariance (c/o-,; — s- J2cr' i^^^''^^'^)aa' ^la' ii be- 
ing the Pauli matrices), and a U(l) orbital symmetry 
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Ci(2)cr.i)- Moreover, we will consider mod- 
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els for which the two legs, or two bands behave identi- 
cally; in other words, we impose a Z2 invariance under 
the permutation of the legs. If we restrict ourselves to 
on-site interactions, the most general model with H = 
U(l)c X SU(2)s X U(l)o X Z2 invariance then reads as 
follows^: 

■i,/cr i i 

i i 

with Hi = Y^icr '^i<y,i {"^laA = 4a.i'^icr,i) being the occupa- 
tion number on the i*^ site. In Eq. ([T|) , the spin operator 
on leg I is defined by 

SL^ = \'^4a.^^af^<^lP,i^ (2) 
a.,13 

whereas = | X^o-C'^io- i ~ '^2<t i) ^^e generator of the 
U(l) symmetry for orbital degrees of freedom. 

Model H]) depends on three microscopic couplings: a 
Coulombic interaction U, a Hund coupling J//, and an 
"orbital crystal field anisotropy" Jt- When Jt ~ 0, the re- 
sulting model is the so-called Hund model which has been 
investigated in the context of orbital degeneracyi^i^i^. 
The generalized Hund model ([T]) is directly linked to ul- 
tracold fermionic -'^''-'^Yb and alkaline-earth atoms with 
nuclear spin / = i^i^. The two-orbital states are de- 
scribed in these systems by the ground state (^Sq = g) 
and a long-lived excited state (■^Pq = e). The almost 
perfect decoupling of the nuclear spin from the electronic 
angular momentum J in the two e, g states ( J states) 
makes the s-wave scattering lengths of the problem inde- 
pendent of the nuclear spin. The low-energy Hamiltonian 
relevant to the ^^^Yb cold gas loaded into a ID optical 
lattice then reads^: 

•Hvb = ~^E (cL,^c/^,i+i + ii.c.) + ^ E "'■^("'/' - 1) 

i.la i.l 
+ y^ng^,ne,^ + Ve:c E 4"/»^e/3,»Cs/3,»Cea,,: > (3) 

where ^ is the fermionic creation operator at site i 
with the nuclear spin-i index a =t,i hi the I = e,g 
electronic states. The occupation number of electronic 
states is nj.i = J2a 4cr i'^icr i- Model ([3]) is then directly 
equivalent to the generalized Hund model ([l} with the 
correspondence: U ^ U + Jt/2, V ^ U - Jt/2 - Jh/4, 
and Vex = -Jh/2. 

Apart from this connection to cold fermions physics, 
one of the main interests of model ([1]) stems from the 
fact that it contains a large variety of relevant mod- 
els with extended continuous symmetries, some of which 
having appeared in different contexts. First of all, when 
Jt = —ZJh/'^, the continuous symmetry is promoted 



to U(l)c X SU(2)s X SU(2)o, and one recovers the so- 
called spin-orbital model^^"— . In absence of Hund cou- 
pling, i.e., Jh = 0, the continuous symmetry group of 
model (HI) is U(l)c x U(l)o x 80(4)^ where each chain 
has a separate SU(2) spin-rotational symmetry. When 
Jj — J///4, model ([l} displays an S0(5) symmetry which 
unifies spin and orbital degrees of freedom. The result- 
ing model, with U(l)c x SO(5)s.o continuous symme- 
try, is relevant to four-component fermionic cold atom 
systems^"—. Finally, when Jh = Jt = 0, one recov- 
ers the U(4) Hubbard model, that has been extensively 
analyzed in recent yearsi^*^"— . As it will be seen in 
Section II, at half-filling, many other highly symmetric 
lines can be identified. For instance, the line Jh = 8U 
unifies spin and charge degrees of freedom with an ex- 
tended U(l)o X S0(5)s,c continuous symmetr}*^. The 
corresponding model has been previously introduced by 
Scalapino, Zhang and Hanke (SZH)^ in connection to the 
SO (5) theory which relates antiferromagnetism to d-wave 
superconductivity^. 

In this paper, we will investigate the nature of the in- 
sulating phases of the zero-temperature phase diagram 
of model ([1]) at half-filling. In this respect, it will be 
shown that the duality approach of Ref. |3l| for half- 
filled fermions with internal global symmetry group H = 
U(l)c X SU(2)^ X U(l)o X Z2 yields eight fully gapped 
phases. These eight Mott-insulating phases fall into two 
different classes. On the one hand, the first class con- 
sists of four doubly degenerate phases which sponta- 
neously break a discrete symmetry of the underlying lat- 
tice model. On the other hand, the second class contains 
four non-degenerate Mott insulating phases. A first one 
is the rung-singlet (RS) phase where two spins on each 
rung lock into a singlet. A second non-degenerate phase 
is a rung-triplet (RT) phase where the spins now com- 
bine into a triplet and this phase is known to be adi- 
abatically connected to the Haldane phase of the spin- 
1 Heisenberg chain^^. Finally, two other Haldane-like 
phases are found: they are spin-singlet but involve two 
different pseudo-spin 1 operators which are respectively 
built from charge and orbital degrees of freedom. In the 
case of the charge pseudo-spin 1 operator, the resulting 
Haldane charge (HC) phase has been found very recently 
in the context of ID half-filled spin-| cold fermions^. 

In addition to this duality approach, it will be shown, 
by means of a one-loop renormalization group (RG) ap- 
proach and numerical simulations (using the Density Ma- 
trix Renormalization Group (DMRG) algorithm^), that 
the zero-temperature phase diagram of model 1^ dis- 
plays seven out of the eight expected insulating phases. 
We find it remarkable that model ([1]), that only has three 
independent coupling constants, turns out to have a rich 
phase diagram which includes the four non-degenerate 
Mott insulating phases. 

Finally, we will make contact with the eight Mott- 
insulating phases found over the years in half-filled gen- 
eralized two-leg ladders with a t± transverse hopping 
terniiiii^. The latter term breaks explicitly the U(l)o 
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symmetry but it is known that this symmetry is recov- 
ered at low-energj*^ii^. The relevant global symmetry 
group is still H and the same duality approach thus ap- 
plies to that case. In this respect, we will connect the two 
families of eight fully gapped phases found for t± = and 
for tj_ 7^ 0. In particular, it will be shown that the two 
problems are in fact connected by an emergent non-local 
duality symmetry. 

The rest of the paper is organized as follows. In Sec- 
tion II, we discuss the symmetries of model ([T]). We also 
present a strong-coupling analysis along special highly- 
symmetric lines which gives some clues about the na- 
ture of the non-degenerate Mott-insulating phases. The 
low-energy investigation is then presented in Section III. 
It contains the duality approach to half-filled fermions 
with internal symmetry group H = U(l)c x SU(2)s x 
U(l)o X Z2. The zero-temperature phase diagram of the 
generalized Hund model ^ and that of highly symmet- 
ric models arc deduced by a one-loop RG analysis. We 
then connect our results to the known insulating phases 
of generalized two-leg ladder models with a t± hopping 
term. In Section IV, we map out the phase diagram of 
model ([1]), SZH, and spin-orbital models with t± = 
by means of DMRG calculations to complement the low- 
energy approach. Our concluding remarks are presented 
in Section V. The paper is supplied with three appendices 
which provide some additional information. Appendix A 
describes the technical details of the continuum limit of 
model The low-energy approach of edge states in 
the non-degenerate Mott-insulating phases are discussed 
in Appendix B. Finally, Appendix C presents the main 
effect of the interchain hopping in the strong-coupling 
regime, close to the orbital symmetric line. 



II. SYMMETRIES AND STRONG COUPLING 

Before investigating the zero-temperature phase dia- 
gram of the generalized Hund model by means of the 
low-energy and DMRG approaches, it is important to 
fully determine the special lines which exhibit enlarged 
symmetry. It turns out that, at half-filling, many highly 
symmetric lines can be highlighted. Their study will give 
some important clues on the possible Mott-insulating 
phases of the model at half-filling. 



A. Highly symmetric lines 

The generalized Hund model ([T]) enjoys a global inter- 
nal symmetry group H = U(l)c x SU(2)s x U(l)o x 
Z2 on top of the lattice discrete symmetries like one-step 
translation invariance, time-reversal symmetry, site and 
link parities. For generic filling, and on four different 
manifolds - that correspond to some fine-tuning of the 
lattice couplings - in the space of coupling constants, this 
model possesses a higher symmetry. 



First of all, on top of the SU(2)s that rotates spin de- 
grees of freedom, one can define a SU(2)o orbital pseudo- 
spin operator: 

Tt = ^K,,-n2,,), (4) 

with ni^i = X^o-"''^.*' I ^ 1,2. When Jt = -3Jh/4:, the 
U(l) orbital symmetry of the Hund model ^ is enlarged 
to SU(2)o, with generators given by Eq. The result- 
ing model displays a U(l)c x SU(2)s x SU(2)o continu- 
ous symmetry and has been considered in systems with 
orbital degeneracy like the spin-orbital mode P^'^^" — . 

A second highly symmetric model is defined for Jh = 
0: then the interacting part of model ([ij simplifies as 
follow: 

+ - y j X!"!'*"'^,^, (5) 

from which we deduce that each leg has a separate SU(2) 
spin rotation symmetry so that the continuous symmetry 
group of model (O is U(l)c x U(l)o x S0(4)s. 

When Jt = JhI^^ as shown in Ref. [SJ, model ([1]) 
is known to be equivalent to the spin-| cold fermionic 
model with interacting part: 

2 

i i m— — 2 

where Ua = {2U - 7Jt)/4, and U2 = {2U + Jt)/4. 
In Eq. (ED, we have P]„^, = Ea/^(^™li, f; a/3)cl,,4,»' 
a, (H — ±|,±i; and (J?n||, |; q;/3) are the Clebsch- 
Gordan coefficients for spin |. Model ([6|) is known to 
exhibit a U(l)c x S0(5)5.o continuous symmetry with- 
out any fine- tuning^. 

Finally, for Jh = Jt = 0, spin and orbital degrees of 
freedom unify to a maximal SU(4) symmetry and model 
(dl) takes the form of the Hubbard model for four com- 
ponent fermions with a U(4) invariance. 

At half-filling, the chemical potential /i is given by 
/io = to ensure particle-hole symmetry. More highly- 
symmetric lines can be found in this particle-hole sym- 
metric case. It stems from the fact that, as in the spin-i 
Hubbard model, the U(l)c charge symmetry can be en- 
larged to an SU(2)c symmetry at half-filling^i^. In this 
respect, one can define a charge pseudo-spin operator by: 

7-t ^ J J _ J J 

Jt - \ {n^~'2), (7) 

which is a SU(2)s spin-singlet that satisfies the SU(2) 
commutation relations. This operator is the generaliza- 
tion in two-leg ladder or two-band systems of the pseudo- 
spin i operator introduced by Andersoi*^ and by Yang 
in eta-pairing problems^. 
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Many interesting lines can then be considered. A sim- 
ple way to reveal them is to consider the energy levels 
of the one-site Hamiltonian ^ with t = Q. The corre- 
sponding spectrum is depicted in Fig. [TJ On top of the 
four symmetric lines that we have identified above, we 
find nine additional lines where higher continuous sym- 
metries emerge (see Table 1). Among all these new highly 
symmetric lines, there are three interesting models with 
two independent coupling constants, i.e., with only one 
fine-tuning. 



Eo= 



ti — T- 



E, = Jt - 2U 



E3= ^-2U 




E4 = 



-3Jh 
4 



2U 



-,4f 

- -H- 



-H- -H- -h ^ 

4f ,- 

- -H- 

-h ^ 



FIG. 1: Energy level diagram for the one-site Hamiltonian ([T]) 
with t = 0. 

A first one corresponds to the SZH model with U(l)o 
X SO (5) continuous symmetry. Such SZH model with 
no transversal hopping t± is defined as follows: Hszh = 



i ^ 

+^szhEK,-1)K,-1) 



As shown in Ref. [a, the fine-tuning Jszh = 4([/szh + 
^SZh) (i-e., Jh = 8U in the context of model ([1])) makes 
the lattice model ([8]) U(l)o x S0(5)s,c symmetric. The 
S0(5) symmetry unifies here spin and charge degrees 
of freedom and is thus different from the spin-| cold 
fermionic atoms ^ case. 

A second symmetric line is found for Jt = 2U with the 
emergence of a SU(2)s x S0(4)c,o continuous symmetry. 
In that case, charge and orbital degrees of freedom play 
a symmetric role and are unified by a SO (4) symmetry. 

Finally, the last extended symmetric ray with the fine- 
tuning Jh ~ — 8J7/3 (see Table I) corresponds to a model 
with U(l)o X SU(2)c X SU(2)s continuous symmetry. 
Such a model has two independent SU(2) symmetries: 
one for the spin degrees of freedom and also a second 
for the charge degrees of freedom. In this respect, it is 
very similar to the spin-orbital model and can be called 
"charge-spin" model. 



B. Strong-coupling analysis 

The identification of these highly symmetric models 
is very useful since several possible insulating phases of 
the generalized Hund model ([1]) can be inferred from a 
strong-coupling analysis. Such an approach has already 
been performed for some special lines of Table I such as 
the half-filled U(4) Hubbard chain^i^, the SO (5) spin-| 
model^i^, and the SZH on(^. 

Here, we present a simple strong-coupling approach 
along three special lines which enables us to identify sev- 
eral non-degenerate Mott-insulating phases. To this end, 
let us first consider the line Jt = 2U with U > and 
Jh < 0. In the absence of hopping term (i.e., t = 0), the 
lowest energy states are the spin triplet (see Fig. [T]). 
An effective Hamiltonian can then be deduced by treating 
the hopping term as a perturbation in the strong coupling 
regime \U, Jh, Jt\ ^ t. At second order of perturbation 
theory, we find an antiferromagnetic SU(2) Heisenberg 
chain: 



He 



Js E (^1 



^2 



^1, 



t+1 



2,i+l 



) , (10) 



Sd, 



where Ca and da are respectively the fermion annihilation 
operator of the upper and lower leg of the ladder with 
spin index a. The occupation numbers on the z*'' site 
are denoted by nc(^d),i respectively. The spin operators 
Sc{d).i a-re defined similarly to those of the Hund model 
(see Eq. ([2])). It is straightforward to relate the SZH 
model to the generalized Hund model ([1]): 

jj _ UsZH + VsZH 

~ 2 
Jh = Jszh 
Jt = UszH — VsZH- 



(8) where Js = -At^/{JH - 4C/) > 0, and Si,i are the spin 
operators defined in Eq. ([2]). The resulting fully gapped 
phase is the well-known RT phase of the two- leg spin-i 
ladder which is adiabatically connected to the Haldane 
phase of the spin-1 chain^^. Such a non-degenerate gapful 
phase displays a hidden antiferromagnetic ordering which 
is revealed by a string-order paramete r^^i^^i^'^ : 

lim {{SI, + S'i,) e''^^'fc='+i(s?,'=+sJ,0 

|i-j|->oo 

x(5f,, + 5y)^0. (11) 

This RT phase is also known to exhibit spin-i edge states 

(9) when open-boundary conditions are considered^^— . 
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TABLE I: Extended continuous symmetries of model ^ at half-filling 



Extended continuous symmetry 


Fine-tuning 


Degenerate levels 


U{l)c X SU{2)o X SU{2)s 


Jt = -3Jh/4: 


i52 = -E4 


U{l)c X U{l)o X SO{4), 


Jh = Q 




U{l)c X SO{5)s,o 


Jt = Jh/4 


-E2 ~ Ez 


U{l)c X SU{4)s,o 




E2 = E3 = E4 


SO (7) 


Jt = 2(7, Jh = 8C/ 


Eq — E2 = E3 


X S[/(4)«,, 


f/ = Jh = 


Eo — Ez = Ei 


517(2), X SO{5)s,o 


Jt = -2(7/3, Jh = -8(7/3 


Eo = E4, E2 = Ez 




Jt — — ULy, JH — OIJ 


Eq = Ez, E2 = E4 


SO(5)c,o X SU{2)s 


Jt = 2(7, Jh = -817/3 


Eo — E2 = E4 


U{l)o X SO(5)s,c 


Jh = 8(7 


Eo = Ez 


SO(4)c,o X SO{4)s 


Jh = 0, Jt= 2(7 


Eo = E2, Ez = Ei 


SO(4),,o X Sf/(2)« 


Jt =2(7 


Eo = i?2 


U{l)o X 5f/(2), X 5f/(2)c 


Jh = -8C//3 


-Eo = E4 



A second interesting line is Jh = —8(7/3 where the 
charge degrees of freedom enjoy an SU(2) symmetry en- 
largement. In the absence of hopping term, the lowest 
energy states for U < and J* > are the Eq , E4 levels 
as it can be seen from Fig. [T] Keeping only these three 
states, we obtain, at second order of perturbation the- 
ory, an effective (pseudo) spin-1 antiferromagnetic SU(2) 
Heisenberg chain: 

"Heft = Jc ^ Jj • Jj-l-i, (12) 

i 

with Jc = -4^^/(6(7- Jf) > 0. The effective Hamiltonian 
(jl2p expresses in terms of the spin-singlet charge opera- 
tor ([7]) which is a pseudo spin-1 operator in the triplet 
states Eq,E4 of Fig. [T] We then expect the emergence 
of Haldanc-like phase for charge degrees of freedom as 
it has been recently found in the context of half-filled 
spin- 1 cold fermions^. Such a HC phase is fully gapped 
and non-degenerate. It displays a hidden ordering that 
is revealed by the string-order parameter: 

lim (Jfe"^i-+i'^^- J,^) 7^ 0. (13) 
\i-j\-yoo 

A deviation from the line Jh = —8(7/3 breaks the 
SU(2) charge symmetry down to U(l) and in the strong- 
coupling regime the lowest correction to model is 
a single-ion anisotropy term DcJ^ii^i)"^ (with Dc = 
3J^f/4 + 2(7). The Haldane phase of the spin-1 chain is 
known to be stable under a weak single-ion anisotropy^. 
A large enough Dc may give rise to a Ising phase with 
(Jf) ^ (i.e. a CDW from the definition 0) or to a 
large- Z?c phase which is a non-degenerate gapped singlet 
phase. The latter, with {J^) = 0, corresponds to the RS 
phase of the two-leg spin- ^ ladder where the two spins of 
the rung bind into a singlet state {E4 state of Fig. [T]) for 
an antiferromagnetic interchain coupling. 

Finally, a last interesting symmetric ray is Jt ~ 
—3Jh/4: where the U(l) orbital symmetry is enlarged 



to SU(2). Along this line, when Jh > and U is not too 
negative, the lowest energy states of the one-site Hamil- 
tonian are levels E2 , E4 of Fig. [T] At second order of 
perturbation theory, we now find a spin-1 antiferromag- 
netic SU(2) Heisenberg chain for the orbital degrees of 
freedom: 

n,s^ JoY,f,-f,+i, (14) 

i 

with Jo = \6t^ /{%Jh + 8(7). We thus expect the emer- 
gence of a new Haldane phase for the orbital degrees of 
freedom that will be called Haldane orbital (HO) phase 
in the rest of the paper. The resulting hidden ordering 
is captured by the following string-order parameter: 

lim (Tfe"^^--+i^'^-T/) 7^ 0. (15) 

\i-]\^oo 

A deviation from the line Jt = —3Jh/4: breaks the 
SU(2) orbital symmetry down to U(l) and in the strong- 
coupling regime the lowest correction to model (fT4l) is 
a single-ion anisotropy term DoJ^ii'^i)'^ (with Do = 
Jt + 3J///4). For sufficiently strong value of \Do\, the 
HO phase will be destabilized into either an orbital den- 
sity wave (ODW) which is described by the E2 states 
with (T/) 7^ or a RS phase, i.e., the E4 state with 
{Tt)=0. 

In summary, the strong-coupling analysis along highly 
symmetric lines reveals the existence of four non- 
degenerate Mott-insulating phases (RT, HC, HO, RS) 
and two gapful phases with long-range density ordering 
(CDW, ODW). 

III. LOW-ENERGY APPROACH 

In this section, we present the details of the low-energy 
approach of the generalized Hund model ([1]) at half- 
filling, in the weak-coupling regime |(7, Jh, Jt\ t. The 
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zero-temperature phase diagram of model (IT} will be in- 
vestigated by means of the combination of a duality ap- 
proach and one- loop RG calculations. In particular, we 
will determine the different insulating phases in the weak- 
coupling regime and make connection with the ones found 
within the strong-coupling approach. 



A. Phenomenological approach 

The starting point of the low-energy approach is 
the linearization around the Fermi points ikp of the 
dispersion relation for non-interacting four-component 
fermions. We thus introduce four left and right mov- 
ing Dirac fermions Li„, Ri„ (Z = 1, 2 and a =t,4.), which 
describe the lattice fermions c/^.j in the continuum limit: 



where the triplet of Majorana fermions ^ = 
accounts for the spin degrees of freedom since the SU(2)s 
spin rotation symmetry of the lattice model ([1]) is gener- 
ated in the continuum by Jir + J2R + L. The Ma- 
jorana fermion is related to the discrete Z2 interchain 
exchange as it can be seen from Eq. (fT8| . Finally, the 
four remaining Majorana fermions can be cast into two 
pairs, each of which is associated to the two U(l) sym- 
metries in H: ^^'^ (respectively ^^'^) Majorana fermions 
describe the orbital (respectively charge) U(l) symmetry. 

With this identification at hand, we can derive the low- 
energy effective theory for the generalized Hund model 
p]) at half-filling. Assuming only four-fermion (marginal) 
interactions, the most general model with H = U(l)c x 
SU(2)s X U(l)o X Z2 invariance can be easily deduced 
from the Majorana fermion formalism: 



i?,<,(x)e*^^" + L,,(.T)e 



— ikpx 



(16) n = -!^^(e^a,e^- ^25.^2) 



with kp = 7r/2ao at half- filling and x = iao (oq being the 
lattice spacing) . The next step of the approach is to use 
the Abelian bosonization of these Dirac fermions to ob- 
tain the low-energy effective Hamiltonian for model ([!]). 
The details of these calculations are given in Appendix 
A. Here, we present a more phenomenological approach 
which is based on the symmetries of the lattice model 

The continuous symmetry of the non-interacting model 
(H]) is SO (8). In the continuum limit, the SO (8) sym- 
metry can be revealed by introducing eight real (Ma- 
jorana) fermions from the four complex Dirac {R^L)i„ 
ones. The non-interacting fixed point is then described 
by the SO (8)1 CFT with central charge c = 4^. The 
chiral currents (1 < a < & < 8), which gener- 

ates this CFT, can be expressed as fermionic bilinears: 
4"^) = ^^L{R)iL(Ry ^here are the eight left (right) 
moving Majorana fermions. 

When interactions are included, the SO (8) symmetry 
is broken down to H = U(l)c x SU(2)s x U(l)o x Z2. 
The key point of the analysis is to identify how the eight 
Majorana fermions of the SO (8)1 CFT act in H. One way 
to obtain the correspondence is to focus on the currents 
which generate the different continuous symmetry groups 
in H. The uniform part of the continuum limit of the spin 
operator ^ on the leg I ^ 1,2 defines the chiral SU(2)i 
currents Jir^l- 

JiL=L]JfL,,,J,R^R\JfR,,. (17) 

As in two-leg spin ladde r^^i^^ , the sum and difference of 
these chiral SU(2)i currents can be locally expressed in 
terms of four Majorana fermions £ ' among the eight 
original ones: 



JlR.L + J2R..L 



S.R,L A ^R^L 



a=7 



■Y.^eRd.CR-eLd.eL) 



IVt 



a=l 
5 



f6 c6 
'iR^L 




+ 



\a=l 

3 

g^Y.^'kCL + giY^^Ril 

a=l a=4 

jZiRCi] +(ekl+ei?i)x 

\a=7 I 

3 5 

9iY.^Rii+ 9^Y.^Rii^ d'^^Ril 

a=l a=4 




(19) 



The diflFerent velocities and the nine coupling constants 
cannot be determined within this phenomenological ap- 
proach based on symmetries. In this respect, a direct 
standard continuum limit procedure of the lattice model 
must be applied. This is done in Appendix A and we find 
the expression of the velocities: 



Vc = Vp + 



ao / 3 



Vs = Vp 



Vt = Vp 



2ti 

Oo 

27r 



U - 
U ~ 



TT \2~' 4 

/ TT Jh_ Jt_ 
2 2 

2 



•^0 (tt , , Jt 



(20) 



JlR,L~J2R,L — ^^R,l£.rL' 



(18) whereas the identification of the nine coupling constants 
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reads as follows 

(tj 3-/^ 
.96 = ao ^3L/ - y ^ 

59-ao(t/-^-|). (21) 

The main advantage of this Majorana fermions descrip- 
tion is that the symmetries of the original lattice model 
are explicit in the low-energy effective model ((T^ in sharp 
contrast to the standard Abelian bosonization represen- 
tation (sec for instance Eqs. (|A5IA7[) of Appendix A 
where the symmetries arc hidden). In particular, using 
Eqs. (pnil2ip . one can check that all extended symmetries 
of Table I are indeed symmetries of model ([TO)) . 

B. Duality approach 

On top of the continuous symmetries of model ([l} , the 
low-energy effective Hamiltonian ([T5)) displays exact hid- 
den discrete symmetries which take the form of duality 
symmetries. Indeed, as shown recently in Ref. [3l|, gen- 
eral weakly-interacting fermionic models with marginal 
interactions exhibit non-perturbativc duality symmetries 
in their low-energy description. Those will help us to list 
and identify possible gapful phases that may occur at 
low-energy. This is the object of the present subsection. 
The duality symmetries are easily identified here within 
the Majorana formalism (|19p since they are built from 
Kramers- Wannier duality symmetries^i^ of the underly- 
ing two-dimensional Ising models associated to the eight 
Majorana fermions they simply take the following 

for m — > — Cli ~^ £,R- Applying the approach of 
Ref. |3l[ yields 8 possible dualities, that can be built out 
of three elementary ones that one may choose as: 

^7 : ei^-ei 

We now describe in detail each of the 8 phases. 



1. Spin-Peierls phase 

The essence of a duality approach is to relate differ- 
ent phases between themselves. In this respect, we need 
a starting phase from which the dual phases can be ob- 
tained. Such a phase can be most simply chosen by con- 
sidering the special fine-tuning VIq : gi = g (i = 1, . . . ,9) 
in Eq. p^ . The RG study presented shortly will assess 
that this line is attractive under the RG flow. On this 
particular line of the nine-dimensional parameter space, 
the interacting part of the low-energy effective model 
takes the form of the SO (8) Gross- Nevcu modcl^i: 

^Pn? = f (^EeSeS^ ■ (22) 

The S0(8) symmetry rotates the eight Majorana 
fermions and is the maximal continuous symmetry of 
the interaction model (HH). The S0(8) GN model is 
integrable and a spectral gap is generated for g > 
Q 70, 72, 73 ^ The resulting phase corresponds to a spin- 
Peierls (SP) ordering with the lattice order parameter 

Csp = EMa(-l)|cL,»C/^,,+i- 

Indeed, a straightforward semiclassical approach to 

model (|22|) reveals that the bosonic fields $c,s./,s/ of the 

basis (|A3p are pinned into the following configurations 

for g > 0: 

($a) = V^Pa, ($a) = ^ + V^^a, (23) 

Paila (a = c,s,f,sf) being integers. In addition, the 
ground state degeneracy of this phase can be deduced 
within this semiclassical approach since there is a gauge 
redundancy in the bosonization procedure (|A1[) : 

*/aflX ^laR,L + \/TTPlaR,L, (24) 

where piaR,L arc integers [l ~ 1,2 and a =tii)- This 
transformation leaves intact the Dirac fermion fields 
Ri„,Li„. Using the change of basis (jASp . we deduce that 
among the field configurations ([23]) . only two of them 
are independent (they cannot be connected by the gauge 
redundancy transformation (|24p): 

($,)=V^, {^sj.sf)=Q. (25) 

These two ground states are related by the one-step 
translation symmetry Toq, which is described in the 
bosonization approach by: $c — >■ ^'c + -v/tt- The SP phase 
is thus two-fold degenerate and spontaneously breaks the 
translation symmetry Tap as it should. 

The continuum bosonized description of the SP order 
parameter is given by: 

Cgp - ]J cos (v^$a) + n ^^"^ (V^*a) , 

a—c,s.f,sf a—c.s,f.sf 

(26) 

from which we deduce that indeed this order parameter 
condenses in the field configurations ([25)1 : (Osp) 7^ 0- 
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2. Charge density wave phase 

Starting from the SP phase, one can infer aU possible 
gapful phases that may appear at low-energy. We ob- 
tain a second degenerate phase hj performing a duality 
transformation fii : — > which is a symmetry 

of model if 37,8,9 —57,8,9- The resulting gapful 
phase, named A^i, is governed by the following interact- 
ing Hamiltonian, which replaces the S0(8) line ([2^ : 



''•int — o 1 / ^ 

\a=l 



'iBSL 



(27) 



In bosonic language, the duality ili only affects the 
charge degrees of freedom and corresponds to a simple 
shift of the charge bosonic field: $cL ^cL + V^/^ and 
^cR ^cR- We deduce from Eq. ([25]) that the phase Mi 
is two-fold degenerate with the two semiclassical ground 
states: 



(28) 



Hence, the Aii phase also breaks the one-step transla- 
tion symmetry. The semiclassical approach enables us to 
identify the Mi phase as a CDW phase, described by 

the lattice order parameter Ocdw = Si,;cr(~l)''^L,i'^;<T,j- 
Indeed, in the bosonization description, the CDW order 
parameter reads as follows: 

CcDW - cos(0F$c) n sin(^/^$a) 

a=sj,sf 

- sin(V^$,) Yl cos(V^$a), (29) 

a=sj,sf 

and it obviously condenses in the ground state configu- 
ration (EHl): (OcDw) 7^ 0. 



3. Orbital density wave phase 



We can define a second duality transformation 
^^'^ — which is indeed a symmetry of model ^ 
if 52.4,8 — ^ —52,4.8- br that case, the S0(8) line 
replaced by: 



'2/O2 



^ I ^ ^rH — Cr^L — ^R^L j 
\a=l,2,3;6,7,8 / 



(30) 



The duality ft2 affects the orbital degrees of freedom 
and is represented by a shift in the orbital bosonic field: 
<i>/L — > ^/L + The resulting phase, M2, is a two- 

fold degenerate with ground state configurations: 

(<!'/) = (*c„s,./)-0, 

(<1>,)=V^, = =0. (31) 



Similarly to the CDW phase, one deduces that M2 has 
a long-range ODW ordering. The order parameter is: 

CoDW = T,^,lai-^yi^^y^^4a,^'^la,i ^hosc bosonized 
form reads: 



CoDW ^ cos(v^$/) Yl sin(V7F$a) 

a—c,s.sf 

- sin(V^$/) Y[ cos(y¥$a), (32) 

a— c,s,s/ 

and it obviously condenses in the ground state configu- 
rations (ED): (OoDw) 7^ 0. 



4. Spin Peierls-n phase 

The last two-fold degenerate phase, M3, which breaks 
translation symmetry, is obtained from the SP phase with 
help of the duality ^3 = ^1^2 : Cl "^'^'® ^ _^4,5,7,8 
It is a symmetry of model (|19p if 52.4,7.9 —52,4.7.9- 
The M3 phase is governed by the following interacting 
Hamiltonian; 



^int - 



2 ( 5Z ^R^L ^ ^r£,L I 

\a=l,2,3;6 a=4,5,7,8 / 



(33) 



The order parameter characterizing the M3 phase is 
a Spin Peierls-TT (SP^r) order parameter with an alter- 
nating dimerization profile on the two legs: Osp^ = 



Its bosonized form is: 



Osp^ ^ n cos(V^$,) Yl sin(V^$,) 

a—c.f a—s,sf 

+ H sin {V^^a) n (\/^*a) - (34) 

a—c,f a—s^sf 

The ground state configurations of the SPjr phase are 
given by: 



(<fc>-^, {^f)^^^ (<!'.,./> =0, (35) 
and (Osp^) 7^ in these configurations. 



5. Haldane charge phase 

So far wc have considered only duality symmetries 
which involve an even number of Majorana fermions. 
A second class of interesting duality symmetries, called 
outer dualities in Ref. [31I is involved when an odd num- 
ber of Majorana fermions is considered. These dual- 
ity symmetries give rise to the second class of Mott- 
insulating phases, the non-degenerate ones, which do not 
spontaneously break the translation symmetry. A first 
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non-dcgencratc gapful phase, named Ai4, is obtained 
from the SO (8) hne 1^ by the duahty : ^^''^'^ ^ 
— Such an outer duahty is a symmetry of model 
(fT9|) when 33,4,7,8 — >■ —33,4,7,8- The resuhing effective 
model for the A^4 phase is: 



SI4 



(36) 



a— 6 



A simple semiclassical analysis of this model shows us 
that the bosonic fields are pinned in the following way: 



i'^c) = V^Pc, ($.,,/) = {Osf) = V^PsJ,sf, 

(37) 

(*c) = ^ + V^Qc, - {Qsf) = ^ + V^1sJ,sf, 

where pa, Qa (a = c, s, /, sf) are again integers. Using the 
gauge redundancy 1^^, we observe that the phase A^4 is 
indeed non-degenerate with ground state configuration: 



(38) 



Unfortunately, the order parameter of this phase 
cannot be written locally in terms of the original lattice 
fermions. In this low-energy procedure, it involves order 
and disorder operators of the underlying two-dimensional 
Ising models. The situation here is similar to the RS and 
RT phases of the two-leg spin-i ladder—. These gap- 
ful phases are non-degenerate and display a hidden an- 
tiferromagnetic ordering and possibly edge states which 
can be revealed through non-local string order parame- 
|-gj. rj26,62^63 ^ In this respect, in order to build those op- 
erators, it is useful to introduce the following quantities 
in terms of the occupation numbers ni^^i of the original 
lattice fermions: 





1 , 


+ ni^.i 


+ n2\.i 


+ n2i.i) , 


^s.z — 


1 1 


- 


+ n2\.i 


- n-2i.i) , 


nf.i = 


1 1 

2 ('^iT,» 


+ ni^.i 


- «2t,i 


- "2;,i) , 


nsf.i = 


1 , 

2 ('^iT,» 


- 


- "2t,i 


+ «2;,i) • 



(39) 

We then consider two classes of string-like order param- 
eters: 



cos 



O 



odd 
a A 




SUa i COS 



(40) 



with Sua^i = Ua^i - {na,i) and a = c,s, f, sf. The string 
operators (|40[) are respectively even or odd under the 



transformation Suaj — >■ —Sua^i- The bosonization de- 
scription of these string-order parameters is cumbersome 
due to the non-locality of the operators in Eq. (|40|) . A 
naive continuum expression can be derived with help 
of its symmetry properties like in the two-leg spin lad- 
one-dimensional extended Bose-Hubbard 

modeli^: 



0° 



dd 



COS (y7r$a) , 



sm 



(41) 



We thus deduce that the odd charge string operator 
displays long-range ordering in the Mi phase: 

hm ) 

- lim (sin(V^$,(x))sin(v^$c(y))) 7^0. (42) 

\x-y\^oo 

Using the charge pseudo-spin operator (O, one im- 
mediately observes that this lattice charge string-order 
parameter is equivalent to the long-range ordering (|13p 
obtained within the strong-coupling approach. We thus 
conclude that the Ai^ phase is a HC phase which is 
adiabatically connected to the HC of the strong-coupling 
approach found in the vicinity of the Jh = — 8C//3 
line. This phase displays a hidden ordering, described by 
Eq. p2)) . and pseudo-spin i edge states, as expected for 
a Haldane phase. Those (holon) edge states carry charge 
but are singlet states as far as the spin and orbital de- 
grees of freedom arc concerned (see Appendix |B]) . This 
result will help us to detect numerically the HC phase in 
the DMRG calculations of Section IV. 



Haldane orbital phase 



The next phase, A^5, is found by applying the dual- 
ity symmetries -^l^'^ with 52,3,8,9 -52,3,8,9- 
The resulting dual interacting Hamiltonian reads 



f25 



E 



Va=l,2,3,7,8 




Its ground state configuration is given by: 



(*/> = 



= (e,s/) = o. 



(43) 



(44) 



The physical properties of the Ai^ phase are very similar 
to the ones for the HC phase. The orbital degrees of 
freedom are central to the Ai^ phase and plays the role 
of the charge degrees of freedom for the HC charge. In 
this respect, the M5 phase is characterized by the long- 
range order of the odd orbital string order parameter: 

lim (01:'^:^0°/f) 



lim (sin 

\x-y\-^OD 



(V^$/(x))sin(y^$/(y))) =^0. (45) 
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Using the orbital pseudo-spin operator ([4]), we find that 
the orbital string-order parameter (|45p is equivalent to 
the long-range ordering (fT5|) obtained within the strong- 
coupling approach. We thus conclude that the A^s phase 
is a HO phase which is adiabatically connected to the 
HO of the strong-coupling approach found in the vicinity 
of the SU(2)o symmetric hue Jt -3Jh/4. The HO 
phase is characterized by pseudo-spin ^ edge states (see 
Appendix |B| which carry orbital quantum number only. 
The orbital edge states will be useful to reveal numeri- 
cally the HO phase by means of the DMRG approach. 



7. Rung-triplet phase 



A new non-degenerate phase, named Ai^, is found by 

12 3 12 3 

applying the duality symmetry fig £,£ ' — > ~S,£ ' to the 
SP phase. The effective interacting Hamiltonian which 
governs the properties of the Aie phase is: 



njfle 9 I \ ^ fa _ \ ^ fa fa 



(46) 



In particular, one observes that no odd string-order pa- 
rameters can condense into the Ai^ phase. The latter 
phase is the standard RS phase of the two-leg spin lad- 
der with antifcrromagnetic interchain coupling that we 
have identified in Section H within the strong-coupling 
approach. This phase has no edge state when open- 
boundary conditions are considered^. 



C. Phase diagram 

The duality symmetry approach, discussed in the pre- 
vious subsection, predicts the emergence of eight insulat- 
ing phases for two-leg electronic ladder with the symme- 
try group H = U(l)c X SU(2)s x U(l)o x Za. However, 
this approach cannot determine which phases actually 
appear in the phase diagram of a particular model like 
P]) with H invariance. To answer this question, we need 
to perform a one-loop RG calculation of model ([T^ with 
initial conditions (|2T|) for the generalized Hund model 



In the A4e phase, the ground state configuration for the 
bosonic fields is: 



{^s) = {esf) = ^, (*c,/)=0, (47) 

from which we deduce that the following string-order pa- 
rameter condenses in this phase: 

hm {O^fO^f) 

~ lim (sin {V^^s{x)) sin (V^$,(y))) ^ . (48) 

\x-y\~^oo 

This order parameter is the standard string-order param- 
eter ([TT]) of the RT phase of the two-leg spin ladder with 
ferromagnetic interchain coupling. The Me phase is thus 
a RT phase with spin-i edge states^. 



8. Rung-singlet phase 

Finally, the last non-degenerate phase, called Mj, is 
obtained from the SO (8) line ([2^ by the duality symme- 
try : ?! ^ -?| 



its effective model is: 



nji^? y I \ ^ fa fa f6 f6 

<a—l,a^6 



(49) 



In the semiclassical approach, the ground state configu- 
ration for the bosons is obtained from the SP one ([25)1 
by changing the bosonic field into its dual Qgf in the 
pinning configuration: 



1. The phases of the generalized Hund model 

The one-loop RG fiow of the nine coupling constants 
of model (|19p can be derived by standard methods^i^. 
By neglecting the velocity anisotropy (i.e. Vc = Vs = 
vq — vt = v) and performing a suitable redefinition of 
the coupling constants (ga — 2t:v fa, a = 1,...,9), we 
find the one-loop RG equations: 



A 




f 2/7^ 




h 


= 2/1/2 + /2/5 + 


hh + 


2/7/8 


h 


= 2/1/3 + 2/2/4- 


h 2/7./9 




h 


= hh + 'if2h + 


2/8/9 




h 


= •ifi + fl + 'ifl 






h 


= 3/I + 2/I + /2 






fr 


= 2/1/7 + 2/2/8- 


h hh - 


^hh 


fs 


= if2h + hh + 


hh + 


hh 


h 




h hh ■ 





($,..,/) -(0,/) =0. 



(50) 



(51) 



The nine coupling constants of the low-energy Hamil- 
tonian are not independent: they depend only on 
the three parameters U, Jh and Jt of the original lat- 
tice Hamiltonian ([T]). We thus need to use the initial 
conditions (pij) to determine which phases do come out 
in the zero-temperature phase diagram of the general- 
ized Hund model. A numerical analysis of these differen- 
tial equations, together with the results of the preceding 
subsection, gives us the phase diagram of the model. We 
find seven insulating phases out of the eight possible ones 
found within the duality approach. The missing phase is 
the SP^ phase. Of course, by adding next-neighbor in- 
teractions to the lattice model ^ without breaking the 
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symmetry group H, the latter phase wih be found. Inter- 
estingly enough, this lattice model with three indepen- 
dent interactions possess the four non-degenerate Mott- 
insulating phases that we have revealed with the help 
the duality symmetry and strong-coupling approaches. 
In Fig. [21 we present a section of the three-dimensional 
phase diagram of the generalized Hund model at U ~ 
—0.005 (t is set to unity) where the seven phases appear. 



0.1 ; 

\ /0.05- 

rt\ ^^cdw 

^^^^ — ^ — L-J^ 




L 

y RS 

\ HC 

1 * 


-0.05 - 

ODW 

-0.1 


o 

X 

- c 

o 
X 

c 


i 0.05 0.1 



CDW-c- ->-HC 




^ ^ U(1)xU(1) 

^ U(1) < > SU(2)2 

FIG. 3: Quantum phase transitions that can occur in the 
generalized Hund model. The letter next to the arrows indi- 
cates which degrees of freedom are critical: c=charge, s=spin, 
o=orbital. 



FIG. 2: Low-energy phase diagram of the generalized Hund 
model at half-filling for U = -0.005t {t = 1). (SP=Spin 
Peierls, CDW=Charge density wave, ODW=Orbital density 
wave, RS=Rung Singlet, HC=Haldane Charge, HO=Haldane 
Orbital, RT=Rung Triplet). At the intersection of the highly 
symmetric lines lie points with even larger symmetry (see Ta- 
ble I). 



([TO]) , describing the physical properties of the spin-| cold 
fermions model, depends on five independent coupling 
constants. The resulting phase diagram, as obtained 
from the one-loop RG calculation, is presented in Fig. 
II 



The duality approach that we used to obtain the phase 
diagram allows for an easy characterization of the quan- 
tum phase transitions. Those transitions are located on 
the self-dual lines, where the coupling constants that 
change their sign when going from one phase to the other 
vanish. Fig. ([3]) summarizes all the phase transitions 
that occur in the phase diagram of the generalized Hund 
model. 

We now present the zero-temperature phase diagram 
of several interesting highly symmetric models with two 
independent coupling constants that we have introduced 
in Section II. 



2. Phase diagram of the S0(5) fermionic cold atoms model 

Let us start with the spin-| cold fermionic atoms 
model ^ which is obtained from the generalized Hund 
model (HI by the fine tuning Jh = iJt- In the continuum 
limit, the coupling constants are naturally parametrized 
by the singlet and quintet pairing Uq and U2' gi = 92 = 
95 = ~Uq - 1/2,93 = .94 = Uo- 31/2,96 = Uo + 51/2,97 = 
9s = 21/2,99 — 2Uo. The effective Majorana model 



f 0.1 y 

RS 

0.05- 

-^°(^ 

-0.1 ^ojo5'~^-^^^^ 


/ y 


-0.05 - 

/ -0.1 


^==^^0^^05 0.1 

CDW 



FIG. 4: Phase diagram of the spin-| cold fermions model at 
half-fiUing (t = 1). 
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ODW xZ/ \\\o.05- 


V SZH 

RS 

S0(7) UszH 


-0.1 -0.05 

SP 

-0.05 - 
-0.1 


\ 

\Sv 0.05 0.1 



coupling constants of model p9|) : 



91 = .93 = -ao ( ^ + y 



2 



92^ 94 = -do U 



95 = -flo U 



96 = do [ 3J7 



3J, 
2 



9^ = 99 = -.92 = -.94 
58 = -gi = -93- 



(53) 



FIG. 5: Phase diagram of the SZH model ([8| with the fine- 
tuning JszH = 4(f/szH + Vszh) at half-filling {t = 1). 



3. S0(5) SZH model 

As mentioned in the introduction, there is a second 
S0(5)-symmetric model embedded in the generalized 
Hund model ([T]): the SZH model In the latter model, 
obtained from ([T]) when Jh = 8U, the S0(5) symmetry 
unifies charge and spin degrees of freedom^. The cou- 
pling constants of the Majorana model for the SZH 
model IHl) read as follows: 



51,6,7 
52,8 
53,9 
94 
95 



2Vs 



ZH) 



OoC/sZH 

-ao (3?7szH + 2Vszh) 
-ao (3[/szH + 4Vszh) 
ao {UszH — 2Vszh) ■ 



(52) 



We obtain the phase diagram shown in Fig. ([5]). The 

cold 
ODW 



phases are very similar to the one in the spin-ij 



fermions model with the substitution: 
and HC ^ HO. 



CDW 



4- Phase diagram of SO (4) models 

We now turn to models which display an extended 
S0(4) symmetry. In Section II, we found two different 
S0(4) models with two independent coupling constants, 
i.e., one fine-tuning with respect to the original general- 
ized Hund model ([1]). When Jh = 0, the lattice model 
JS]) enjoys aU(l)c x U(l)o x S0(4)s continuous symme- 
try. The phase diagram of this model can be determined 
by the low-energy approach from the identification of the 



The resulting phase diagram is presented in Fig. |6l 
It contains an interesting line Jt = 2U with enlarged 
S0(4)c,o X S0(4)s symmetry (see Section |TT|. It is 
straightforward to see that along this line for [/ < 0, 
where the transition between the CDW and ODW ap- 
pears, the spin degrees of freedom are gapped while the 
charge and orbital degrees of freedom are critical. Hence, 
the quantum phase transition is described by a S0(4)i 
CFT with central charge c = 2. 



o.ih 
0.051 

CDW 1 

1 1 


- / SP 

/ SU(4)3„xU(1)<, u 
^====- \ ^ 


-0.1 -0.05 / 
0/ 

^/-0.05- 

A/ 

/ 

/ -0.1 


^ _ao5 0.1 

X 

J' ODW 



FIG. 6: Phase diagram of the S0(4) model ^ with Jh = 
{t — 1) at half-filling. The continuous symmetry group is 
U(l)c X U(l)o X S0(4),. 



The second model is defined for Jt = 2U (see Table I) 
with S0(4)c.o X SU(2)s continuous symmetry. Here, the 
SO (4) symmetry unifies the charge and orbital degrees of 
freedom. The coupling constants of the continuum limit 
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of this model are given by 



of the continuum hmit of this model read as follows: 



gi ^ -ao \ 2U - 
Jh 

.92 = 97 = flo- 
ss = -ao (2C/ + 



■In 

2 



Jh 



'SJh 

.94 = .99 = ^'^0^— 
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(54) 



from which we deduce the phase diagram of Fig. [71 The 
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FIG. 7: Phase diagram of the half-filled S0(4) model with 
Jt = 2U {t = 1) which unifies charge and orbital degrees of 
freedom. 



quantum phase transition between RS and RT phases for 
Jh = and [/ > is now governed by the spin degrees 
freedom and the S0(4)i CFT. 



The phase diagram of the spin-orbital model is presented 
in Fig. [51 In particular, it includes two non-degenerate 
Mott-insulating phases: HO and RT phases. 
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FIG. 8: Phase diagram of the spin-orbital model at half-filling 
with Jt = -3Jif/4 {t = 1). 



Phase diagram of spin- orbital (charge) models 



The generalized Hund model (jlj reduces to the spin- 
orbital model when Jt = — 3Jh/4 with a U(l)c x SU(2)o 
X SU(2)s symmetry which has been studied in the con- 
text of orbital degeneracy^"—. The coupling constants 



Finally, we consider the related spin-charge model 
which is defined for the fine-tuning Jh = — 8J7/3 with 
U(l)o X SU(2)c X SU(2)s continuous symmetry at half- 
filling. The resulting model is similar to the spin-orbital 
model where now spin and charge are put on the same 
footing. We find the following continuum limit for the 
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spin-charge model: 
91 
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(56) 



Its phase diagram is depicted in Fig. [HI Here, two non- 
degenerate Mott insulating phases appear: the HC and 
RT phases. 
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FIG. 9: Phase diagram of the spin-charge model at half-filling 
with Jh = -SU/3 [t = 1). Fhe symmetry is: SU{2)s x 
S'C/(2)c X U{l)o. 



D. Effect of an interchain hopping term 

Let us now consider the generalized Hund model ([T]) 
to which we add an interchain hopping term: 



'la. i^2cr, i 



h.C, 



(57) 



This problem has been previously studied and its phases 
are known^iiiii^. We will show that all those phases 
(there are eight of them) can be connected one-by-one 
(by a non-trivial duality) to the phases of the general- 
ized Hund model ([ij . Our approach will use a weak cou- 
pling analysis, and our conclusion relics on the existence 



of a hidden orbital symmetry emerging at low-energy. In 
Appendix \C\ we show that this approach is also valid 
in the strong coupling regime, close to the orbital line 
Jt = -3Jh/^- 

In the presence of interchain hopping, the non- 
interacting part of the Hamiltonian cannot be readily 
diagonalized. We need to do a change of basis and intro- 
duce bonding and antibonding operators: 



JO- 



i(ci. + (-i)^c2,). 



V2 



(58) 



Using this basis, the kinetic part can be diagonalized in 
momentum space; there are now two decoupled bands 
(bonding and antibonding bands) and two Fermi points 
(provided tj_ < 2t) kpi and kp^ such that, at half-filling, 
kpi + kp^ = Tr/ao. One could proceed by linearizing 
around the four Fermi points, by introducing continu- 
ous bonding and antibonding fermionic fields Raa{x) and 
Laa{x) (with a = 1, 2, (7 =t,i) by bosonizing, refermion- 
izing, and expanding the interactions in this basis. In- 
stead, we will follow an approach that is based on the 
(continuous) symmetry content of the theory, by showing 
that the symmetry of the low energy, continuous theory 
is essentially the same for ij^ 7^ as for t±^ = 0. 

This is not obvious: naively, the interchain hopping 
term breaks the U(l)o orbital symmetry and the anal- 
ysis of the preceding sections breaks down. However, 
as noticed before^ii^, at low energy another U(l) sym- 
metry emerges in the orbital sector. To see this, let us 
consider the difference S = kpi — kp2 = 2 arcsin(t_L/2i). 
It is a continuous function of t±, and provided d is not 
commensurate to tt, in the continuum limit and at weak 
coupling, one can safely ignore umklapp terms that os- 
cillate at wave-vector that are integer multiples of S. 
Retaining only marginal, four-fermions interactions of 

the form Hi^i^Qi^o-i, with VP^*^ = R,L, one sees that 
in order for this term to be non-oscillating and to give 
a contribution to the interacting part of the continu- 
ous theory, it has to conserve separately the quanti- 
ties p+ = NiR + N2L and p_ = Nil + N2R, where 
NaL = Y.aS dx Ll^L^^ and NaR = EaIdx RLRacr are 
the total number of left (right) fcrmions in each of the 
bonding and antibonding bands. It results that the differ- 
ence p+ — p_ = {NiR — N2r) — {Nil — N2l) is conserved: 
this is nothing but (twice) the total orbital current in the 
bonding/antibonding basis along direction z. This quan- 
tity generates a U(l) orbital symmetry 11(1)^. One thus 
concludes that the low-energy continuous theory has a 
symmetry U(l)cXSU(2)s x 11(1)^ x Z2, the same has in 
the tj^ = case. 

The total orbital current — p- can be mapped onto 
the total orbital charge (all other conserved quantities, 
i.e., the electric charge and the SU(2)s spin generators, 
being unaffected) by the duality U,±_ that leaves all Majo- 
rana fcrmions {£,°')a^b invariant, but changes ^£ to — ■Ci- 
The duality Q.± is highly non-local. A little algebra shows 
that it has the following action on the bonding and anti- 
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bonding modes 



1,2): 



"at 



L 



t 



L 



a 4, 



-L 



t 



(59) 



With those elements at hand, the general form of the 
low-energy Hamiltonian can thus be readily deduced from 



([T9| by performing the replacement — >■ —(X- While 
this approach tells nothing about the (bare) value of the 
coupling constants ga , it tells us that the structure of the 
RG equations (|5ip is the same, and is thus sufficient to (i) 
enumerate the phases of the model and (m) relate each 
of those phases to the phases of the generalized Hund 
model ([U. 

From the preceding sections (i_L =0), we deduce that 
there are eight possible insulating phases. The effect of 
the duality ^± is to transform the bosonic field into 

its dual with a shift of ^/2: 6/ + ^/2. We 

then obtain the pinnings of the bosonic fields from those 
found for t_L = via the duality fij,; they are summa- 
rized in Table HIl In the language of Ref. [HI VL± is an 
outer duality (it affects only one Majorana fermion, or 
equivalently, it maps a bosonic field onto its dual): it re- 
sults that it maps degenerate phases onto non-degenerate 
ones and vice-versa. 

The pinnings of Table |TT] allow us to identify each phase 
with one of the already known phases of the generalized 
two-leg Hubbard laddei— The SP phase becomes an 
S-Mott phase with order parameter: 



O 



S-Mott 



(60) 



The CDW phase becomes the D-Mott phase, which is 
the same phase as RS, and is characterized by the order 
parameter: 



O 



D-Mott — ^ Ci^,iC2l,i 



(61) 



The ODW (respectively SP^r) phase give an S'-Mott (re- 
spectively D'-Mott) state which differs from the S-Mott 
(respectively D-Mott) only in the pinning of the charge 
bosonic field $c (see Table |lT|^. Such order parame- 
ters have the slowest decaying correlation function when 
the system is doped. The non-degenerate RS phase will 
become the two-fold degenerate CDWtt phase with an 
inter-leg phase differenced^: 



CDW-ir 



i,l(T 



(-ir(-i) 



'+ic^ c 

1(7. i 1(7 A ' 



(62) 



The HC phase becomes a p-density wave (PDW) phase 
which is described by the condensation of the order pa- 
rameter—: 



OPDW = ^(-i)x-i)'+' (c; 



(63) 



The HO phase gives a staggered flux (SF) phase (or a d- 
density wave phase) whose ground states display currents 



circulating around a plaquctte, with order parameter: 



(64) 



This phase spontaneously breaks time-reversal symme- 
try. Finally, the RT phase will become a f-density wave 
(FDW). This phase consists in currents flowing along the 
diagonals of plaquettes: 



(65) 



It also breaks time-reversal symmetry. 

As already emphasized, the emergent duality symme- 
try VL±^ provides only a correspondence between the set of 
phases with = with those oit±^ 7^ 0. By no means, 
it maps a given model deflned by {t±^ = 0, U, Jh, Jt) onto 
the model (tj_ 7^ 0, U, Jh, Jt)- In this respect, it is nec- 
essary to use the one-loop RG calculation to map out 
the phase diagram of the generalized Hund chain with a 
transverse hopping. Using Eq. ([51]) and the initial con- 
ditions for that model 
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(66) 



we find the presence of the eight insulating phases. Fig- 
ure [in] presents a section of the phase diagram for U = 
— 0.005i where the eight phases are revealed. Interest- 
ingly enough, the generalized Hund chain model with a 
transverse hopping turns out to be the minimal model 
for two-leg electronic ladder models which contains the 
eight insulating phases found within the low-energy ap- 
proacbiiii^. In particular, one does not need to add fur- 
ther interactions (first neighbors density interactions for 
instance) contrarily to Refs. [T1IIT2I . 
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TABLE II: Pattern of the bosonic fields pinning for t± — and t± 7^ 0. 
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FIG. 10: The eiglit fully gapped phases of the generalized 
Hund chain with a transverse hopping at half-filling {U — 
-0.005t). 



IV. DMRG CALCULATIONS 

We now carry out numerical calculations using DMRG 
in order to investigate the various phase diagrams. For 
simplicity, we will restrict ourselves to the t± ~ case 
for which wc can use a purely one-dimensional imple- 
mentation of the model. Moreover, for a system of fi- 
nite size L, we can fix the total number of particles 
N = 2Qc = 2 ''^c,i = 2L, the 2-component of the total 
spin 5*^ = J2i ''^s,i as well as the z-component of the to- 
tal orbital operator = J2i''^f,i (^^^ Eq. ([5^ ). so that 
the states are labeled by the triplet {Qc, jT''). Typi- 
cally, we keep up to 1600 states, which allow to have an 
error below 10~^, and we use open boundary conditions 
(OBC). 

As it has been revealed by the low-energy approach, 
there are eight possible insulating phases, that are re- 
lated by duality transformations. In order to character- 
ize them, we can either compute local quantities (bond 
kinetic energy, local density, . . . ) to identify states that 
break a given symmetry, or investigate the presence or 
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FIG. 11: Numerical phase diagram of the generalized Hund 
model at half-filling for IJ — —t. Notations are the same as 
in Fig. [3 



absence of various edge states to detect non-degenerate 
Mott phases. Thus, the three possible Haldane phases 
that have been predicted (RT, HC, and HO) can be 
characterized by looking for the presence of edge states 
with quantum numbers (Qc Sz,Tz) respectively equal to 
(L, 1,0), {L + 1,0,0) and (L, 0,1) (see details in Ap- 
pendix |B|) . 

From these measurements, wc can draw various cuts 
of the phase diagram in the (?7, Jh, Jt) parameter space. 
In Fig. [n] we present our data for fixed U = —t. Seven 
(out of eight) insulating phases are found and the overall 
topology nicely agrees with the low-energy predictions 
shown in Fig. [51 although those were obtained with a 
much weaker interaction U = — 0.005t. As another ex- 
ample, we consider a highly symmetric case, namely the 
SZH model (see Section IIII C 3|) . which corresponds to 
the case Jh = 8U. We recall that the low-energy phase 
diagram of this model is shown in Fig. \5\ The numeri- 
cal phase diagram (Fig. [T2|) is in excellent agreement not 
only for the overall topology but also for the location of 
the phase boundaries. 

Finally, we also did simulations for the spin-orbital 
model obtained by fixing Jt = —iJu/^- In that case 
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FIG. 12: Numerical phase diagram of the SZH S0(5) model 
with the fine-tuning Jszh = 4:{Uszh + Vszh) at half-filling. 
Notations are the same as in Fig. [S] Dashed lines indicate 
models with higher symmetry (see Table 1). 



as well, we observe that the low-energy predictions (see 
Fig. |S]) and our numerical data (see Fig. [T31) coincide 
very well. 



Jfj 



-4 



. ♦ 



▲ 
▲ 



♦/ 

4 



* 

_i_ 



▲ 


CDW 






RT 




♦ 


SP 




► 


HO 






1 




2 


3 





u 



FIG. 13: Numerical phase diagram of the spin-orbital model 
(Jt — — 3Ji^/4) at half-filling. Notations are the same as in 
Fig. |8l Dashed lines indicate models with higher symmetry 
(see Table 1). 

In brief, for all the cases we considered when t± = 0, 
wc obtain an excellent agreement between the low-energy 
predictions (obtained for very weak couplings) and the 
numerical phase diagram (obtained at moderate cou- 
plings). This is an important result of our paper. 

The case of a transverse hopping t± ^ is much more 
difficult to analyze with respect to the low-energy ap- 
proach. As discussed in Section III D, in the latter ap- 
proach, there is an emergent U(l) symmetry in the or- 
bital space which is not present on the lattice. It means 
that we need to consider long-size systems in the DMRG 



calculations in order to compare with the low-energy pre- 
dictions. A second difficulty arises with respect to the nu- 
merical discrimination between the S-Mott (respectively 
D-Mott) phase and the S'-Mott (respectively D'-Mott) 
phase. A numerical analysis of relevant string-order pa- 
rameters is clearly called for to distinguish them. To the 
best of our knowledge, wc are not aware of any numeri- 
cal study which reports the existence of the (S,D)'-Mott 
phases in generalized two-leg electronic ladders. We plan 
to investigate elsewhere the numerical phase diagram of 
the generalized Hund model with a transverse hopping. 



V. CONCLUDING REMARKS 

In summary, we have investigated the phase diagram 
of the generalized Hund model ([T]) with global symmetry 
group H = U(l)c X SU(2), x U(l)o x Za at half-ffiling 
by means of complementary low-energy and DMRG tech- 
niques. The interest of this model is manifold. First, it 
covers and unifies several highly symmetric lattice mod- 
els: the S0(5)-symmetric model describing spin-| cold 
fermions; another (different) S0(5)-symmetric model in- 
troduced by Scalapino, Zhang and Hanke that unifies 
antiferromagnetism and d-wave superconductivity; and 
two SU(2)xSU(2)-symmetric spin-orbital models. While 
its orbital U(l) symmetry can seem a little odd from the 
point of view of applications to the description of Hub- 
bard two-leg ladders, we show that at weak-coupling, it 
shares the same continuous low-energy effective theory 
with the well known Hubbard two-leg ladder with in- 
terchain hopping. Third, it is directly relevant to the 
description of Ytterbium 171 loaded into an optical ID 
trap. We are able to treat on an equal footing all the 
phases appearing in those models coming from differ- 
ent contexts. The rich phase diagram of the generalized 
Hund model has to be contrasted with the minimal char- 
acter and simplicity of this model, which only depends 
on three microscopic couplings. 

Wc bricfiy recall our main results: by means of a du- 
ality approach, wc predict that the phase diagram for 
half-filled four-component fermions with global symme- 
try H = U(l)c X SU(2)s X U(l)o X Z2 consists in eight 
Mott-insulating phases. These phases fall into two differ- 
ent classes. A first class consists of two-fold degenerate 
fully gapped density phases which spontaneously break 
a discrete symmetry present on the lattice. The sec- 
ond class comprises four non-degenerate Mott-insulating 
phases which arc characterized by non-local string or- 
der parameters: RS, RT, HC, and HO phases. A one- 
loop RG calculation for model ([1]) reveals the existence 
of seven phases out of the eight ones consistent with the 
duality approach. The missing phase, which can be found 
by adding further nearest-neighbor interactions, is an al- 
ternating bond ordered phase (SPtt). These results have 
been confirmed numerically by means of a DMRG ap- 
proach for moderate couplings. In this respect, we found 
an excellent agreement between the two complementary 
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approaches. 

Finally, we have connected our low-energy results to 
the eight previously known insulating phases found in 
generalized two-leg ladders with a transverse hopping tj_ 
term. When t± ^ 0, the U(l) orbital symmetry is lost 
on the lattice but becomes an emergent symmetry at 
low-energy. The duality approach with global symmetry 
group H can then still be applied in the presence of a in- 
terchain hopping tj_ term. In this respect, we discovered 
a non-local duality which maps the eight Mott-insulating 
phases for t± = onto the eight phases previously known 
for tj_ ^ 0. The one-loop RG approach to the generalized 
Hund model with interchain hopping predicts the stabi- 
lization of the eight Mott-insulating phases. The latter 
model is thus the minimal model for two-leg electronic 
ladders which displays the eight Mott-insulating phases 
at weak coupling. We also discuss the fate of this emer- 
gent U(l) orbital symmetry when going from the weak- 
coupling to the strong coupling regime. We show that the 
weak- and strong-coupling phases coincide in the vicinity 
of the orbital symmetric line, but our analysis suggests a 
breakdown of this picture away from this special line. 

We hope that future experiments on ^''^Yb or alkaline- 
earth cold fermions atoms loaded into a ID optical lattice 
will reveal some of the exotic insulating phases found in 
our study. 



tion relation: 



[<^laR, $/'<t'l] = -^Sw6„„, . (A2) 

The presence of the Klein factors Kia ensures the correct 
anticommutation of the fermionic operators. The Klein 
factors satisfy the anticommutation rule {k/o-, k/zo-'} = 
2<5/;'^CTcr' and they are constrained so that = 1, with 
F = Ki-\Ki^K2-]-K2i- Hereafter, we will work within the 
F = 1 sector. It will be convenient to work with a pair 
of dual non-chiral bosonic fields: ^jo- = ^i^l + ^laR and 
0/cr = $/crL — ^laR- Last, let us introduce a SU(4) basis 
that will allow us to separate charge and non-Abelian 
(spin) degrees of freedom: 
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In sharp contrast with incommensurate fillings, there is 
no spin charge separation at half-filling. Indeed, in this 
special case, chiral umklapp processes couple those de- 
grees of freedom. Consequently, the resulting low-energy 
Hamiltonian corresponding to model ([T]) takes the form: 

^ = 'He + T~Ls + Humklapp- 

The charge degrees of freedom are described by: 



Appendix A: Continuum Limit 



In this Appendix, we present the technical details of 
the continuum limit of the generalized Hund model ([l} . 
The starting point of the low-energy approach is the lin- 
earization around the Fermi points -iikp of the dispersion 
relation for non-interacting fermions. Four left and right 
moving Dirac fermions Li„, Ri„ {I = 1,2 and a =^, \.) are 
then introduced to describe the lattice fermions Q^.i in 
the continuum limit (see Eq. 

The next step of the approach is the introduction of 
four chiral bosonic fields ^iaR,L through the Abelian 
bosonization of Dirac fermiona^i^: 



Kla 



exp iV 47r$;c,H 



cxp (-i\/47r<I>(o-L ), 



(Al) 



where vp = a^t is the Fermi velocity. 

The part of the bosonized Hamiltonian corresponding 
to the non-abelian degrees of freedom is: 



a.=sj,sf 
a=sj,sf 



cy^'^c//" -f As cy*"c 



5^/ 



(A5) 



where the bosonic fields satisfy the following commuta- 



where we used the compact notation: Cf = cos(/3$a). 



19 



and = cos (/30a)- The bare parameters are: 



Majorana fermions through: 
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Finally, the umklapp part of the Hamiltonian reads: 



^ AeC^cf^ + ATCf^C^ 

+ A^cf^C^ + AgCf^G^ , (A7) 



with 
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(A8) 



In the end, a refermionization procedure will allow us 
to make the exact U(l)c x SU(2)s x U(l)o x Z2 contin- 
uous symmetry explicit in the effective Hamiltonian. For 
this purpose, we introduce eight left and right moving 
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(A9) 



where 771,2,3,4 are again Klein factors which ensure the 
adequate anticommutation rules for the fermions. Using 
this correspondence rules, equations (|A41 IA5IA7P can be 
expressed in terms of these eight Majorana fermions. We 
thus finally obtain the low-energy effective theory for the 
generalized Hund model dD) at half-filling: 
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where the different velocities and couplings are given by: 

Jh , Jt 



on the Dirac fermionic fields of Eq. (1161) are thus 



"0 ( TT 



ao ( jj 3Jt 



51 = -ao ( [/ - — + y 

52 = -ao ([/- — - y 
-an\U + — + — 



3Jh Jt 
~ "2 

2 

2 



97 = ao{U+ — - — 



54 = -ao yU + 
95 = -ao (^U - 
56 = ao I 3J7 - 



59 = ao U7 



3Jh Jt 



(All) 



Appendix B: Edge states 

In this Appendix, we investigate the possible existence 
of edge states in the non-degenerate gapful phases Hal- 
dane charge, orbital and spin, when open boundary con- 
ditions are used. The nature of the edge states will rein- 
force our characterization of, respectively, the HC, HO, 
and RT phases as being pseudo-spin 1 chains in, respec- 
tively, the charge, orbital, and spin degrees of freedom. 
Since the occurrence of edge states in the RT phase has 
already been discussed at length in Ref. [H, we restrict 
our attention in the following to the HC and HO phases. 



1. Open boundary formalism 

The OBC are taken into account by introducing two 
fictitious sites and 4- 1 in Eq. ([T]) and by imposing 
vanishing boundary conditions on the fermion operators: 
Co = Cat+i = O^^Ti^. The resulting boundary conditions 



Lia (0) 



-Rla (0) 



L) 



(Bl) 



with L ^ {N + l)ao and ^ = 1, 2, cr =t,|. The left and 
right-moving Dirac fermions are no longer independent 
due to the presence of these open boundaries. Using the 
bosonization formula (|Aip . we deduce the boundary con- 
ditions on the chiral bosonic fields: 



<^iaL (x = L) 



-'^laR. {X = L) 



L 

ao 



(B2) 



where pio- and gjo- are integers. The total bosonic field 
with internal degrees of freedom la thus obeys Dirichlet 
boundary conditions: 



^i^ {x = L) 



L 

ao 



(B3) 



The next step of the approach is to introduce the mode 
decomposition of the bosonic field compatible with 
these boundary conditions: 



$/<T {x,t) = 



sin (mrx/L) 

n=l 



mr 



X 

L 











2 \ao J 





vp-t/ L 



h.c. 



(B4) 



where Tioia is the zero-mode operator with spectrum 
s/iiqio and anio is the boson annihilation operator obey- 
ing: [a„;o-, a]^;,^,] = 6n,mSia,i'a' ■ The mode decomposi- 
tion of the dual field &ia can then be obtained from the 
property: dtOia = vpdx^ia- 



Qia {x,t) = $o;<T 
cos (niTx/L) 



^ f L 



2 yao 

imrvpt/L 



- 2 + Uola 



h.c. 



(B5) 



with [I>o;ct, Horcr'] = i&ia,i'o'- In particular, and 
= dxQiiy satisfy the equal-time canonical com- 
mutation relation for bosons: [^(^(i, a;), H/zo-' (<, y)] = 
iSia-^i'a-'SLix — y), 5l{x) being the delta function at fi- 
nite size: 5l{x) = e"'^^/^/2L. Using the defini- 
tions ^iaB.,L ~ {^icr ± 'dia)/2, the mode decomposition 
of the chiral bosonic fields ^iaR,L can then be deduced 
from Eqs. (|B4IB5I) . One can then show that these chiral 
fields satisfy the following commutation relations when 
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['^laR.hix) ,'^l'a'R.L{y)] ^ T ^(5,^,/ V Sgll (o: - ?/) (B6) 

1 

sgn(x) being the sign function. At this point, one should 
note a technical subtlety which will play its role in the 
investigation of the possible edge states of the Haldane 
phases. When considering OBC, the sign of the commu- 
tator between [^iaR,^i'a'L] is the opposite of the bulk 
one (|A2[) . The latter comes from the identity often used 
in the bosonization approach^: 

•^laR.L (x) = i {x) tJ dy n,, (y)^ , (B7) 

which does not take properly into account the boundary 
conditions on the fields. This subtlety has no effect on 
the derivation of the low-energy Hamiltonian (|19p which 
is still valid in presence of OBC as it can be easily shown. 
However, it will be important for the discussion of edge 
states as first observed in Ref. [6^ for the determination 
of boundary excitations of the semi-infinite two-leg spin 
ladder. 

With this formalism at hands, we are now in posi- 
tion to investigate the possible edge states in the gen- 
eralized Hund model ([1]) with OBC. To simplify the dis- 
cussion, we will consider a semi-infinite geometry where 
the OBC is located at the i = site and L — )■ -|-oo. 
The low-energy effective Hamiltonian density is still given 
by Eq. ([T^ . but now we see, using the bosonic bound- 
ary conditions (jB2[) . the change of basis (jASp . and the 
refermionization (jA9p . that the eight Majorana fermions 
^R L = Ij ■ • • 1 8) must verify the following boundary 
conditions: 



(0) 



(0). 



(B8) 



2. Edge states in the Haldane charge phase 

Let us set our investigation on the symmetry line 
Jh = — 8C//3 where the charge degrees of freedom display 
an extended SU{2) symmetry and form the pseudo-spin 
1 in charge. For strong attractive U, we expect that 
the spin and orbital gaps will be larger than the charge 
gap. This is confirmed numerically by DMRG calcula- 
tions: for instance, for U/t ^ —3 and Jt = 0, the spin 
and orbital gaps are close respectively to 6t and 5t, while 
the charge gap is roughly O.lt. Within this hypothesis, 
we can safely integrate out spin and orbital degrees of 
freedom of model p9)) . For simplicity, let us choose to 
work on the line Jt ~ Jh/4:, where the spin and orbital 
degrees of freedom are unified. Keeping only the rele- 
vant terms and again neglecting velocity anisotropy, the 



resulting effective Hamiltonian then reads: 

8 

c 

a— U " ^ 

8 

^ f a ^ 



n 



im. 



dx iad.a - ClO^Cl) 
/ dx J2CrCl, (B9) 

•^0 a=6 



with the mass TOc given by (57 ~ on the line Jh 
-8C//3): 



W7 



(BIO) 



Model (|B9p is the sum of three decoupled semi-infinite 
free massive Majorana fermion models. Hence, let us now 
consider a single Majorana fermion ^r^l model: 



'^Lrn = --Z- I dx {^Rda;£.R- £.Ldx^L)-'im / dx^R^L, 

^ Jo Jo 

(Bll) 

with boundary condition: ^l(O) = £,r{0). Model (IBllP is 
quadratic with dispersion relation ek = V v^k? + and 
with the fermionic decomposition^: 

£.R{x,t) = -^J2{^k{cos{kx + 0k) + ism{kx))e-'^''' 



k>0 



h.c.} 



9 (m) e 



£,L{x,t) = -7= {Cfc (cos {kx + 9k) ~ i sin {kx)) e" 

V 2/7 , ^ ^ 



fe>0 



h.c.} + 



6 (to) e' 



-rnx j V 



(B12) 



where ^k is fermion annihilation operator with wave- 
vector fc, 9 is the step function, and 77 is a zero mode real 
fermion normalized according to 7]"^ = i. In Eq. (jB12P . 
9k is defined by: 



exp {i9k) = 



vk - 



(B13) 



The key point of Eq. ()B12p is the existence of an expo- 
nentially localized Majorana state with zero energy inside 
the gap (midgap state) for a positive mass to. In con- 
trast, for negative to, such a zero mode contribution does 
not occur since it is not a normalizablc solution. 

The presence of edge states for the HC phase thus de- 
pends on the sign of the mass rric. Using the definition 
(IA9P and the commutator (IB6P. we find: 



TOc 



97 
TTflo 



2 \^-f 



(B14) 



In the HC phase, we have: {^s) = ($/) = (6s/) = 
(see Table II), so that TOc = —5gr/2TT. Hence, since 
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gr = aoU, the mass rric > 0, which signals the emergence 
of three locahzed Majorana modes 77" (a = 6, 7,8) from 
the mode decomposition ()B12p . Moreover, three local 
Majorana fermion modes are known to define to a local 
pseudo spin-i operator J"- thanks to the identitj^: 

^" = -^6°^+^^+', (B15) 

that is a consequence of the anticommutation relations 
{77", rj^} = (5"''. Wc thus conclude on the existence, in the 
HC phase, of a pseudo spin-i edge state at the boundary 
which can be viewed as a holon edge state. 

One recognizes that the pseudo-spin projection along 
= —irfrf is proportional to the total charge 
generating U(l)c: in the continuum limit, within 
the convention ()A9|) . one has^ Qc = j^i^i — > 

/ dx {e^en + iiii) = - 2z E.>o {ae: - ae:) , 

showing that the zero-mode contributes the total charge. 
In a finite size system of size L with two boundaries, the 
edge states come into pairs, that organize into a pseudo- 
spin singlet and a pseudo-spin triplet. Edge states at the 
two end of the chain interact, leading to a singlet /triplet 
splitting that goes to zero in the thermodynamical limit. 
It results that one observes a mid-gap state with quan- 
tum numbers {Qc, S\T^) = (L ± 1, 0, 0). 

3. Edge states in the Haldane orbital phase 



the charge and spin degrees of freedom: 

TTflo 

+ \{C^ + C^))- (B18) 

so that TOc = —Sgs/STT in the HO phase (see Table II). 
As on the considered line, — —'2aoU, the mass nio > 
0. We thus conclude on the existence of three localized 
Majorana modes which form a pseudo spin-i (orbital) 
edge state at the boundary in the HO phase. 

Repeating the same argument as in the HC phase leads 
to the expectation - in a finite geometry with size L and 
two OBC - of a mid-gap state with quantum numbers 
(Qc,5^^-) = (L,0,±l). 

Appendix C: Strong coupling around the orbital line 

In this Appendix we give a description of the effect of 
an interchain hopping tj_ in the strong coupling regime. 
We will first show that close to the orbital symmetric 
line {Jt = —'3Jh/4: ), a low-energy effective continuous 
theory can be derived and trivially solved, allowing for 
the identification of the phases and the phase transitions. 

Then, wc will show that a U(l) symmetry emerges at low 
energy, similarly to what is known at weak coupling (see 
Section III D). 



Let us now sit on the symmetry line Jt = — 3Jh/4 
where it is the orbital degrees of freedom that display 
an extended SU(2) symmetry and form a pseudo-spin 1. 
For the sake of simplicity, let us choose to look at line 
Jh = 8U when charge and spin degrees of freedom are 
unified into an SO (5) symmetry. For repulsive U, we ex- 
pect that the charge and spin gap will be higher than the 
orbital gap (which is confirmed numerically by DMRG 
simulations for instance for U/t = 1: the charge and spin 
gaps are both equal to 6. it while the orbital gap is O.OSi) 
so that we can safely integrate out the corresponding de- 
grees of freedom. The resulting leading effective Hamil- 
tonian is (neglecting velocity anisotropy): 




with the mass ruo given by: 

5 

mo = ^9sY.(^'hCL)■ (B17) 
0=1 

The latter mass can be expressed in terms of the 
ground state expectation value of the bosonic fields for 



1. Continuous theory 

The effect of the inter-leg hopping in the strong cou- 
pling regime is in general a complicated problem: since 
the hopping term breaks the U(l) orbital symmetry, this 
process will induce transitions amongst on-site states (de- 
picted in Fig. [ij that belong to different symmetry mul- 
tiplets. 

However, as noticed in Section fll Bl for a special fine- 
tuning of the couplings Jt = — 3Jh/4 the lattice model 
enjoys a SU(2) symmetry in the orbital sector. Close to 
this line, the orbital symmetric line, a strong coupling ex- 
pansion can be performed: orbital degrees of freedom are 
the only low-energy modes, and an effective Hamiltonian 
for the orbital operators can be derived, that governs 
their dynamics. Noticing that the interchain hoping term 
can be expressed in terms of orbital degrees of freedom, 
t± Ea,»(cL,«C2^,i + cL,iCi^,i) = 2ij. E» Tf, one remarks 
that the effect of t± close to the SU(2)o line amounts to 
the analog of a transverse magnetic field in direction x 
for orbital degrees of freedom, resulting in the following 
effective Hamiltonian: 

Heff = JoY.T,- f,+i + DoY^iT^-f + kY, t:, (ci) 

i i i 

where = 16tV(9JH + 8U), Do = Jt + ^Jn/i. and 
ho = 2t±. Orbital operators Ti being spin-one operators. 
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close to the orbital line the problem is thus equivalent 
to a spin-1 Heisenberg model with single- ion anisotropy 
under a transverse magnetic field. 

In the absence of a magnetic field ho, it is known that 
the spin-1 Heisenberg chain with a single- ion anisotropy 
can be described in terms of continuous degrees of free- 
dom, namely 3 Majorana fermions^^ (to be consistent 
with the main text, we call them a = 4, 5, 6), that are 
related as follows to the uniform components of the spin 

operators: T^" ^ T'^_(x) = + ^^^''^'^r") + 

. . ., where the dots indicate oscillating terms. The low 
energy spectrum of the theory is well described, at low- 
est order, by a theory of 3 free massive fermions, corre- 
sponding to 3 branches of magnons. Now the magnetic 
field also yields a term that is quadratic in fermions. Ne- 
glecting velocity anisotropies, we thus end up with the 
following quadratic Hamiltonian: 



-He 



+ i m„CS^£ + tho {iUl + iUV) -(02) 



The masses arc phenomenological parameters. When 
Do ~ t± = 0, one has a single mass scale m oc 
Jo (the gap of the spin-1 Heisenberg chain), and in 
the general case one can parametrize them as 7714 = 
7715 = m — S and me = m + with S (x Do at 
first order in Do/ Jo- One immediately sees that the 
fermion decouples. Fourier transforming the remain- 
ing fermions, ^2(i?)(^) ~ I intro- 
ducing *fc ^ (ei(fc),e|(fc),ei(fc),C^(fc)) the Hamilto- 
nian reads X)fc>o ^ -kUk^k- The one-particle spectrum 
is obtained from the eigenvalues of the matrix: 



Uk = 



/ 


Vok 




iho 


\ 






-Vok 





iho 




-iho 





Vok 


irriQ 


V 





-iho 




-Vok J 



(C3) 



This yields two branches with energies: e± (k) 

mPS'^, where eo(fc) 



el{k) + hl + S'^ ± 2yjhlel{k) + m^d^, where eo[lt) is 

the spectrum at Do = t± — 0: eo{k) — ^Jv^}i^^^\^Tr? . 

We can now identify several phase transition lines, for 
which the spectrum is massless, i.e., admits modes of ar- 
bitrary low energy. A first Ising transition line is readily 
obtained when S = m: then = 0. 

Other critical lines are found by solving the equation 
e_(fc*) = (the branch e+(fc) is always gapful): one finds 
that for ho+S^ = m^, this equation has a solution k* ~ 0, 
and for 5 = and \ho\ > m, two solutions k* = ±|fc* | 7^ 0, 



with vik* 



On the former critical line, the 



massless degrees of freedom consist in a single Majorana 
mode and one has a central charge c = | . The latter crit- 
ical line has central charge c — 1, and corresponds to the 
well-known commensurate/incommensurate transition of 



the isotropic spin-1 Heisenberg chain under a magnetic 
field. 



5/m 



RS=D-Mott 



HO] \ 



SF 



h/m 





FIG. 14: (Color online) Phase diagram of the generalized 
Hund ladder with interchain hopping. The coordinate on 
the horizontal axis is proportional to t± whereas that on the 
vertical axis is the distance Do to the SU(2) orbital sym- 
metric model. Blue lines indicate c = i critical lines, and 
the red one is a c = 1 critical line. See the main text for 
the definition of the phases. The dashed line, with equation 
m^S^ — ha{hi — m^), indicates the location where incommen- 
suration appears: on the right of this curve, the lower band 
e_ displays a minimum for a wave vector kmin- The two in- 
sets display the typical spectrum in the commensurate and 
incommensurate regions, respectively. 

To identify the phases of the ladder lying on both sides 
of those lines (see Fig. [T4| . one notices that by continu- 
ity with what happens at t_L = 0, the phase with 7714 < 
(obtained at large enough 5) must be a RS phase (which 
coincides with the D-Mott phase). The phase at 7714 > 
is readily identified by means of the duality transforma- 
tion — > that amounts to $/ o — 0/ (see 



Eqs. dSSl)): it is a SF phase. 

By continuity with t± = 0, the pocket (5^ 4- /i^ < m? 
corresponds to a HO phase. The last phase to be mapped 
out is the one lying at -\- 5'^ > w? and 5 < 0: it is 
adiabatically connected to the ODW phase that one has 
obtained at tj^ = and Do < so that it is an ODW 
or CDW^ phase, as we named it in the context of the 
ladder with an interchain hopping. 



2. Effective low-energy theory 

We now make a qualitative connection with the weak- 
coupling analysis of Section IIII Dl which crucially relies 

on the existence of an extended U{1) symmetry in the 
orbital sector, and show that this symmetry is also emer- 
gent at low energy in the strong-coupling regime. In the 
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strong-coupling limit, in the vicinity of the orbital sym- 
metric line, the interchain hopping t± amounts to a mag- 
netic field in orbital space, that will eventually drive the 
system to a incommensurate state: when (5 = 0, as soon 
as ho > ni the Fermi points shift to kp = + k* . 
Now if the "magnetic field" is large enough, i.e., if k* is 
large enough, this leaves the space for a low-energy de- 
scription built on fields expanded around the new Fermi 
points. In this effective low-energy description, Akp umk- 
lapp terms have no effect as being strongly oscillating. 
This remains true if one departs from the orbital sym- 
metric line (5 = 0: in this case, the spectrum of the low 
energy band, e-(fc), develops a gap but the incommen- 
suration is still present in the form of a minimum (which 
becomes infinitely deep in the limit ho ^ m) in e_ (fc) at 

a value 't^o^mm — 

,/hl — 771,2(1 + j^). This happens (see 



Fig. [14]) when 6 is smaU enough, \S\ < 5c = ho\l ^ — 1- 
But is there an associated emergent symmetry, as it was 
the case in the weak-coupling limit? 

To investigate this, one has to enter into further de- 
tails. Denoting by the unitary matrix diagonalizing 
Uk, with AlUkAk = Diag(e_(A:), -e_(fc), e+(fc), -e+(fc)), 

and introducing ~ — (x- (A), X+(fc)), one 

can represent in the cigenmodc basis the orbital U(l) 
generators along .t, Tf^^ ^ = i J dx {^Ul±^U%) = 



Ik>o ^k^k T^ot^± are the total Noethcr charge and 

current. The Hermitian matrix elements L arc too 
cumbersome, and not particularly enlightening, to ap- 
pear here, but the matrix •* has a remarkable block 
structure: 



{-)^( h{k)I t2{k)a^ +h{k)a^ 

\t2{k)cj' +h{k)cj'' -h{k)l 



Let us now introduce two quantities characterizing the 
spectrum: the absolute minimum of the lower band, 
and the gap A from the minimum of the lower band Eq 
to the upper band (see insets of Fig lT4| . Now, if the 
two bands e+(fc) and e_(fc) are well separated, i.e. if 
A ^ £^01 it makes sense to describe the theory at low 



energies in terms of the two Majorana fermions X-(^)- 
Introducing the projector P~ that projects on this sub- 
space, the effective Hamiltonian iJ^ = P^HP^ = 

(fc)X- (^)'''zX- (fc) commutes with P~ T^ot - P - 
Thus, the total orbital current along x is asymptotically 
conserved in the limit of large A/i?o j and the low-energy 
theory is effectively U(l) symmetric. For a physical quan- 
tity computed in the low-energy, U(l)-symmetric theory, 
violation of this emergent symmetry by processes con- 
necting the two bands will result in corrections of order 
(i?o/A)2. 

Now, the parameters Eq and A bear qualitatively dis- 
tinct forms in the commensurate and incommensurate 
regions. For simplicity, we restrict our attention to the 
case ho > m. In the commensurate region, one has: 



Eq 

A 



\/ + hi — m, 



2m, 



(C5) 



so that in general Eq/A is not small. 

In contrast, in the incommensurate region, one has: 



A 




(C6) 



so that the symmetry U(l) is violated by terms of order 

{6/2ho)'^: the symmetry U(l) is asymptotically exact in 
the limit ho ^ S. One thus concludes that close to the 
line Jt = — 3Jh/4, the situation at large coupling coin- 
cides with that at small coupling. One can thus reason- 
ably infer that this symmetry is emergent in all regimes at 
least close to the line Jt = — 3Jff/4. Our strong coupling 

analysis suggests that the emergent U(l) is broken when 
moving far away from the orbital symmetric line ; how- 
ever, the question of the existence of such an emergent 
symmetry in the strong-coupling regime of the generic 
electronic two-leg ladder goes far beyond the scope of 
this paper. 
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